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Abstract. The classic middle-thirds Cantor set leads to a singular continuous measure via a 
distribution function that is know as the Devil's staircase. The support of the Cantor measure 
is a set of zero Lebesgue measure. Here, we discuss a class of singular continuous measures 
that emerge in mathematical diffraction theory and lead to somewhat similar distribution 
functions, yet with significant differences. Various properties of these measures are derived. 
In particular, these measures have supports of full Lebesgue measure and possess strictly 
increasing distribution functions. In this sense, they mark the opposite end of what is 
possible for singular continuous measures. 
■ For each member of the family, the underlying dynamical system possesses a topological 

factor with maximal pure point spectrum, and a close relation to a solenoid, which is the 
Kronecker factor of the system. The inflation action on the continuous hull is sufficiently 
explicit to permit the calculation of the corresponding dynamical zeta functions. This is 
achieved as a corollary of analysing the Anderson-Putnam complex for the determination of 
the cohomological invariants of the corresponding tiling spaces. 



Dedicated to Robert V. Moody on the occasion of his 70th birthday 

1. Introduction 

The probably most widely known singular continuous measure emerges as the unique in- 
variant probability measure for the iterated function system [31J of the classic middle-thirds 
Cantor set. The construction and the distribution function F of the resulting measure are 
illustrated in Figure [TJ Due to its shape, F is known as the Devil's staircase. It is a non- 
decreasing continuous function that is constant almost everywhere, which corresponds to the 
fact that the underlying measure gives no weight to single points, but is concentrated on an 
uncountable set of zero Lebesgue measure (the Cantor set). The Cantor measure is thus both 
continuous and singular, hence purely singular continuous. 

Singular continuous measures occur in a wide range of physical problems, most notably 
in the theory of non-periodic Schrodinger operators; see |20[ [23] and references therein for 
examples. In particular, it is an amazing result that singular continuous spectra are in a 
certain sense even generic here; compare [HJ [39]. One would also expect the appearance 
of singular continuous measures in mathematical diffraction theory |22l [30] [8], where the 
Thue- Morse sequence provides one of the few really explicit examples. Recent experimental 
evidence |54] indicates that this spectral type might indeed be more relevant to diffraction 
than presumed so far. This case has not yet received the theoretical attention it deserves, 
though partial results exist in the dynamical systems literature; compare [35] 132], 
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Figure 1. The distribution function F of the classic middle-thirds Cantor 
measure. The construction of the Cantor set is sketched in the inset. 

The Thue- Morse system is an example of a bijective substitution of constant length [45J. 
This class has a natural generalisation to higher dimensions, and is studied in some detail 
in |26} 127] : see also |55| and references therein for related numerical studies. Bijective sub- 
stitutions form an important case within the larger class of lattice substitution systems. For 
the latter, from the point of view of diffraction theory, a big step forward was achieved in 
|36 } 137 1 [38]. where known criteria for pure pointedness in one dimension [23] were generalised 
to the case of Z d -action. Moreover, a systematic connection with model sets was established 
(see [311 132] f° r detailed expositions and [35] for a rather complete picture), and there are also 
explicit algorithms to handle such cases; compare [281 E] an d references therein. Nevertheless, 
relatively little has been done for the case without any coincidence in the sense of Dekking 
[23] or its generalisation to lattice substitution systems [361 EZ]- Although it is believed that 
one should typically expect singular continuous measures for bijective substitutions without 
coincidence, explicit examples are rare. 

As a first step to improve this situation, we investigate a class of generalised Thue-Morse 
sequences in the spirit of [34]. They are defined by primitive substitution rules and provide 
a two-parameter family of systems with purely singular continuous diffraction. Below, we 
formulate a rigorous approach that is constructive and follows the line of ideas that was 
originally used by Wiener [53] , Mahler [40] and Kakutani [32] for the treatment of the classic 
Thue-Morse case. Some of the measures were studied before (mainly by scaling arguments 
and numerical methods) in the context of dimension theory for correlation measures; compare 
[331 ES] and references therein. 
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The paper is organised as follows. We begin with a brief summary of the Thue-Morse 
sequence with its spectral and topological properties, where we also introduce our notation. 
Section [3] treats the family of generalised Thue-Morse sequences from [6], where the singular 
continuous nature of the diffraction spectra is proved and the corresponding distribution 
functions are derived. Here, we also briefly discuss the connection with a generalisation of 
the period doubling sequence. The latter has pure point spectrum, and is a topological factor 
of the generalised Thue-Morse sequence. This factor has maximal pure point spectrum. The 
diffraction measure of the generalised Thue-Morse system is analysed in detail in Section [H 
and its Riesz product structure is derived. In Section [5l we construct the continuous hulls of 
the generalised Thue-Morse and period doubling sequences as inverse limits of the substitution 
acting on the Anderson-Putnam cell complex [3], and employ this construction to compute 
and relate their Cech cohomologies. The substitution action on the Cech cohomology is then 
used in Section [6] to derive the dynamical zeta functions of the corresponding substitution 
dynamical systems. Finally, we conclude with some further observations and open problems. 

2. A SUMMARY OF THE CLASSIC THUE-MORSE SEQUENCE 

The classic Thue-Morse (or Prouhet-Thue-Morse, abbreviated as TM below) sequence [2] 
can be defined via a fixed point of the primitive substitution 

(1) Q = Qtm : 1^11, 1^11 

on the binary alphabet {1,1}. The one-sided fixed point starting with 1 reads 

(2) v = v v x v 2 ... = 11111111 ... , 

while v is the fixed point starting with 1. One can now define a two-sided sequence w by 

Vi, i > 0, 
v-i-i, i < 0. 



(3) Wi 



It is easy to check that w = . . . w_ 2 w_ 1 \w w_ 1 ... = ... v 1 v Q \v v 1 . . . defines a 2-cycle under 
g, and hence a fixed point for g 2 , with the central seed 1|1 being legal. Recall that a finite 
word is called legal when it occurs in the n-fold substitution of a single letter for some n € N. 
An iteration of g 2 applied to this seed converges to w in the product topology, which is thus 
a two-sided fixed point of g 2 in the proper sense. 

The sequence w defines a dynamical system (under the action of the group Z) as follows. 
Its compact space is the (discrete) hull, obtained as the closure of the Z-orbit of w, 

X™ = {ShL \i € Z}, 

where S denotes the shift operator (with (Sw)i = Wj+i) and where the closure is taken in the 
local (or product) topology. Here, two sequences are close when they agree on a large segment 
around the origin (marked by |). Now, (X™, Z) is a strictly ergodic dynamical system (hence 
uniquely ergodic and minimal [451 152] ) . Its unique invariant probability measure is given via 
the (absolute) frequencies of finite words (or patches) as the measures of the corresponding 
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cylinder sets, where the latter then generate the (Borel) cr-algebra. Its minimality follows 
from the repetitivity of the fixed point word w, which also implies that X™ = LI(w), which 
is the local indistinguishability class of w. The latter consists of all elements of {1,1} Z that 
are locally indistinguishable from w. 

Here, we are interested, for a given w € X™, in the diffraction of the (signed) Dirac comb 

(4) u = w5 z := w n 5 n , 

where the symbols 1 and 1 are interpreted as weights 1 and —1. This defines a mapping from 
X™ into the signed translation bounded measures on Z (or on R). Since this mapping is a 
homeomorphism between X™ and its image, we use both points of view in parallel without 
further mentioning. 

Given any w € X™, the autocorrelation measure of the corresponding u exists as a conse- 
quence of unique ergodicity. It is defined as the volume-averaged (or Eberlein) convolution 



7 = u) ® oj = lim 



N^oo 2N+1 

where u N is the restriction of u to [-N, N] and Ji is the 'flipped-over' version of the measure 
fj, defined by Jt{g) := fJ,(g) for continuous functions g of compact support, with g(x) = g(—x). 
We use this general formulation to allow for complex weights later on. A short calculation 
shows that the autocorrelation is of the form 

with the autocorrelation coefficients 

N 



(5) ^)= 7V 1 ^^vTiE 



., , - , / , w n W n — m . 

N^oo 2N+ 1 ^ 
n=-N 

Note that 7 applies to all sequences of X™, by an application of the ergodic theorem, as 
77(771) is the orbit average of the continuous function w 1— > uiQW_ m and our dynamical system 
is uniquely ergodic. 

Remark 1 (Alternative approach). Without the diffraction context, it is possible to directly 
define the function rj : Z — > C by (0). It is then a positive definite function on Z, with a 
representation as the (inverse) Fourier transform of a positive measure [i on the unit circle, by 
the Herglotz-Bochner theorem [27]. Which formulation one uses is largely a matter of taste. 
We follow the route via the embedding as a measure on M, so that we get 7 = \i * 5 Z together 
with its interpretation as the diffraction measure of the Dirac comb uj. 

We can now employ the special structure of our fixed point w to analyse 7. One finds that 
r/(0) = 1, r](—m) = 77(771) for all m € Z and 

N-l 



(6) 77(m)=lim-X 



Af->oo TV 

n=0 
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for all m > 0. Here, the structure of w and its relation to v was used to derive from ([5]). 
Observing that v satisfies v 2n = v n and v 2n+ i = v n for all n > 0, one can employ © to infer 
the linear recursion relations 

(7) r)(2m) = rj(m) and rj(2m+l) = — — (j](m) + 77(771 +1)) , 

which actually hold for all m G Z. These well-known relations [32] will also follow from our 
more general results in Section [3] as a special case. One finds r/(±l) = —1/3 from solving the 
recursion for m = and m = —1 with r/(0) = 1, while all other values are then recursively 
determined. 

To analyse the diffraction measure 7 of the TM sequence (following [40|. [32] ) , one can start 
with its pure point part. Defining U(N) = EjL-jvfaC")) > one derives £(4iV) < §£(2iV) 
from the recursion ([7]); see [5] for the detailed estimate needed. This implies jjU(N) jV ~ s> °°> 0. 
By Wiener's criterion [53], this means (7) = 0, so that 7 is a continuous measure (see 
Wiener's Corollary in [33] Sec. 1.7.13] or Wiener's Lemma in [43 \ Sec. 4.16] for details). 

Defining the (continuous and non-decreasing) distribution function F via F(x) = j([0, x]), 
another consequence of ([7]) is the pair of functional relations 




Splitting F into its sc and ac parts (which are unique and must both satisfy these relations) 
now implies backwards that the recursion ([7J) holds separately for the two sets of autocor- 
relation coefficients, r/ sc and r/ ac , with yet unknown initial conditions at 0. Since this means 
Wadl) = ~^Vac(fy together with r] ac (2m) = rj ac (m) for all m G N, an application of the 
Riemann-Lebesgue lemma forces ?7 ac (0) = 0, and hence rj 3C (m) = for all m G Z, so that 
also (7) = 0; compare |32j . This shows that 7 is a singular measure. With the previous 
argument, since 7 7^ 0, we see that it is a purely singular continuous measure. Figure [2] shows 
an image, where we have used the uniformly converging Volterra-type iteration 

1 f 2x 

K+Ax) = 2 J ( x ~ cos M)K{y) d y with F o( x ) = x 

to calculate F with sufficient precision (note that F(x + 1) = F(x) + 1, so that displaying F 
on [0, 1] suffices). In contrast to the Devil's staircase, the TM function is strictly increasing, 
which means that there is no plateau (which would indicate a gap in the support of 7); see 
[S] and references therein for details. 

Despite the above result, the TM sequence is closely related to the period doubling sequence, 
via the (continuous) block map 

(8) <f>: 11, 11 1-^0, 11, Hi-)- 6, 

which defines an exact 2-to-l surjection from the hull X™ to X pd , where the latter is the 
hull of the period doubling substitution defined by 

(9) £> pd : ai-^ab, b^aa. 
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Figure 2. The strictly increasing distribution function of the classic, singular 
continuous TM measure on [0, 1]. 

Viewed as topological dynamical systems, this means that (X pd ,Z) is a factor of (X™,Z). 
Since both are strictly ergodic, this extends to the corresponding measure theoretic dynamical 
systems, as well as to the suspensions into continuous dynamical systems under the action 
of R. The dynamical spectrum of the TM system then contains Z[^] as its pure point part 
|45| , which is the entire dynamical spectrum of the period doubling system. We thus are in 
the nice situation that a topological factor with maximal pure point spectrum exists which is 
itself a substitution system. 

The period doubling sequence can be described as a regular model set with a 2-adic internal 
space [TBI ES] and is thus pure point diffractive. As another consequence, there is an almost 
everywhere 1-to-l mapping [50J [IS] of the continuous hull (see below) onto a (dyadic) solenoid 
S = §2- Here, a solenoid S m (with 2 < m £ N say) is the inverse limit of the unit circle under 
the iterated multiplication by m. The dyadic solenoid is obtained for m = 2. 

The discrete hull X™ of the TM sequence has a continuous counterpart (its suspension), 
which we call Y™. Instead of symbolic TM sequences, one considers the corresponding 
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tilings of the real line, with labelled tiles of unit length. Such tilings are not bound to have 
their vertices at integer positions, and the full translation group K acts continuously on them. 
The continuous hull of a TM tiling is then the closure of its M-orbit with respect to the local 
topology. Here, two tilings are close if, possibly after a small translation, they agree on a 
large interval around the origin. For the same reasons as in the case of the discrete hull, the 
corresponding topological dynamical system (Y™,R) is minimal and uniquely ergodic, so 
that every TM tiling defines the same continuous hull. Similarly, a continuous hull is defined 
for the period doubling sequence. 

According to [3], the continuous hull of a primitive substitution tiling can be constructed as 
the inverse limit of an iterated map on a finite CW-complex T, called the Anderson-Putnam 
(AP) complex. The cells of T are the tiles in the tiling, possibly labelled to distinguish different 
neighbourhoods of a tile within the tiling. Each point in T actually represents a cylinder set 
of tilings, with a specific neighbourhood at the origin. The substitution map on the tiling 
induces a continuous cellular map of T onto itself, whose inverse limit is homeomorphic to 
the continuous hull of the tiling. 

Let T n be the CW-complex of the nth approximation step. It is a nice feature of the 
corresponding inverse limit space Y = hni T n that its Cech cohomology H k (Y) can be com- 
puted as the direct (inductive) limit of the cohomologies of the corresponding approximant 
spaces H k (T n ). In our case, we have a single approximant space F, and a single map (the 
substitution map g) acting on it. Consequently, H k (Y) is the inductive limit of the induced 
map g* on H k (T). Analo gous inverse limit constructions also exist for the hull of the period 
doubling tiling, Y pd , and for the dyadic solenoid, S. 

As a consequence of the above, there is a 2-to-l cover cf>: Y™ — > Y pd , and a surjection 
ifi: Y pd — > S which is 1-to-l almost everywhere. These maps induce well-defined cellular maps 
on the associated AP complexes; see also [19] for a general exposition. We represent these 
maps by the same symbols, <f> and ijj. They induce homomorphisms on the cohomologies of 
the AP complexes, so that we have the following commutative diagram: 

ff fc(r«d) H k (TP d ) H k (T™) 

(10) x2 e >* e* 

#fc( r soi) H k (TP d ) H k (r™) 

All these maps are explicitly known. The inductive limits along the columns not only give 
the cohomologies of the continuous hulls, but also determine the embeddings under the maps 
<jf and -0*. Although ff 1 (Y pd ) and i7 1 (Y™) are isomorphic, the former embeds (under (/>*) 
in the latter as a subgroup of index 2, which reflects Y™ being a two-fold cover of Y pd . 
Furthermore, we get iJ 1 (Y pd )/0*(ii" 1 (S)) = z - % an application of pEJ Prop. 4], compare 
also \18\ Ex. 7] j this corresponds to the fact that there are exactly two orbits on which the 
map ip fails to be 1-to-l. These two orbits are merged into a single orbit under ip. 

The action of the substitution on the cohomology of the AP complex H k (T), more precisely 
the eigenvalues of this action, can be used [3] to calculate the dynamical zeta function of 
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the substitution dynamical system, thus establishing a connection between the action of the 
substitution on H k (T) and the number of periodic orbits of the substitution in the continuous 
hull. We skip further details at this point because they will appear later as a special case of 
our two-parameter family, which we discuss next. 

3. A FAMILY OF GENERALISED THUE-MORSE SEQUENCES 

The TM sequence is sometimes considered as a rather special and possibly rare exam- 
ple, which is misleading. In fact, there are many systems with similar properties. Let us 
demonstrate this point by following [341 [6]. where certain generalisations were introduced. In 
particular, we consider the generalised Thue-Morse (gTM) sequences defined by the primitive 
substitutions 

(11) Q = Q k ,f 1^1^, 

for arbitrary k,£ € N. Here, the one-sided fixed point starting with v = 1 satisfies 

I v m , if < r < k, 



(12) v, 



m(k+l)+r 



if k<r < k + l, 



for m > 0, as can easily be verified from the fixed point property. A two-sided gTM sequence 
w can be constructed as above in Eq. ([3]). Also analogous is the construction of the discrete 
hull X fc | as the orbit closure of w under the Z-action of the shift, where X 1 1 = X™. We will 
drop the index when it is clear from the context. 

Proposition 1. Consider the substitution rule g= g kl of Eq. (fTTI) for arbitrary, but fixed 
k,£ E N. The bi-infinite sequence w that is constructed from the one-sided sequence v of 
(fT2|) via reflection as in Eq. ([3]) is a fixed point of g 2 , as is w' = g(w). Both sequences are 
repetitive and locally indistinguishable. 

Moreover, w is an infinite palindrome and defines the discrete hull X = {S l w \ i S Z} under 
the action of the shift S. This hull is reflection symmetric and minimal, and defines a strictly 
ergodic dynamical system. Similarly, when turning w into a tiling of R by two intervals of 
length 1 that are distinguished by colour, the closure Y of the M-orbit of this tiling in the local 
topology defines a dynamical system (Y, R) that is strictly ergodic for the M-action. 

Proof. Most claims are standard consequences of the theory of substitution dynamical systems 
[45] . as applied to g. The fixed point property includes the relation g 2 (w) = w together with 
the legality of the central seed w_ 1 \w around the marker. Note that each fixed point of g n , 
with arbitrary n E N, is repetitive and defines the same hull X. The latter is minimal due to 
repetitivity, and consists precisely of the LI class of w. Since w and w' coincide on the right of 
the marker, but differ on the left in every position, neither can have a non-trivial period (this 
would contradict their local indistinguishability) . This is an easy instance of the existence of 
distinct proximal elements |17] , Consequently, X = LI(u>) cannot contain any element with a 
non-trivial period, so that w and hence g is aperiodic. 
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The action of the shift is clearly continuous in the product topology. Unique ergodicity 
follows from the uniform existence of all pattern frequencies (or from linear repetitivity). 
This means that (X, Z) defines a topological dynamical system that leads to a strictly ergodic 
dynamical system (X,Z, u), where the unique measure v is defined via the frequencies of 
patches as the measures of the corresponding cylinder sets. The claim about the extension 
to the M-action on Y follows from the suspension of the discrete system, which is easy here 
because the constant length of the substitution g implies that the canonically attached tiling 
is the one described. □ 



Let us mention in passing that the discrete hull X is a Cantor set, while the local structure 
of the continuous hull Y is a product of an interval with a Cantor set; compare [20] and 
references therein. 

Since each choice of k,£ leads to a strictly ergodic dynamical system, we know that all 
autocorrelation coefficients (as defined by Eq. ([5]), with 1 = — 1) exist. Clearly, we always 
have 7/(0) = 1, while several possibilities exist to calculate f?(±l) = fifef- 

As before, we turn a gTM sequence w = (wn) n ^x t ne Dirac comb 

(13) u) = ^2w n 5 n , 

which is a translation bounded measure. Its autocorrelation is of the form 



(14) 7 = 5 m 

rneZ 

with the coefficients rj(m), which can alternatively be calculated via the one-sided fixed point 
v as in Eq. ([6]). Let us now derive a recursion for r](m). Since this will be the 'golden key' for 
almost all spectral properties, we provide a detailed proof. 

Lemma 1. Consider the gTM sequence defined by the primitive substitution g of (|lip . for 
fixed k,£ £ N. When realised as the Dirac comb of Equation (|13p . each element of the 
corresponding hull X = X fc ^ possesses the unique autocorrelation 7 of (|14p . where the auto- 
correlation coefficients satisfy n(0) = 1 and the linear recursion 

r]((k+e)m + r) = — -j(a kAr r?(m) + a kAk+e _ r r/(m+l)) , 

with a k ir = k + £ — r — 2 mm(k, £,r,k + £ — r), valid for all m G Z and < r < k + £. In 
particular, one has n((k+£)m) = rj(m) for m € Z. 

Proof. Existence and uniqueness of 7 are a consequence of Proposition [H via an application 
of the ergodic theorem. The support of the positive definite measure 7 is obviously contained 
in Z, so that 7 = r]5 z . 
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Table 1. Conditions on r and t for the splitting of the sum in Eq. (|15p . and 
the resulting term for the sum from the recursion relation (|12p . 



Conditions 


Term 


< r < k 


< t < k-r 
h — t <C t <" mini A 1 h-\-f — v) 
min(A;, k+i — r) < t < k 


Vri V n -\- m 
^n-\-in-\-l 


k < t < max(fc, k+i—r) 
max(fc, k+i-r) <t< k+l 


V n V n +m+l 


k <r < k+i 


< t < min(£;, k+l-r) 
min(A;, k+i— r) < t < k 


Vn "Vn+m+1 


k < t < max(fc, k+l— r) 
max(fc, k+i-r) <t< 2k+i-r 
2k+i-r < t < k+i 


Vn v n+m+l 
V n V n +m+l 



Since r/(0) = 1 is immediate from Eq. ([5]), it remains to derive the recursion. We begin 
with m > and use formula (El). When r = 0, one finds 



N-l j k+i-lN'-l 

j((k+i)m) =lun QW J2 v nVn+(k+ i )m =^(^^7 E E ' 

_ k+t-lN'-l ( . v „ . , , ^ 

©r 1 \ - \ - I w n u n+m , if < t < k I 

" iV™oo (k + i)N' ^ yv n v n+m , if k < t < k+ij 



k+i-l N'-l 

iV'Too (A; + i)N : 



linl E E = 7 ?( 771 )' 



t=0 n =0 

where the penultimate step follows because v n v n +m = v n v n +m due to E {±1}- 
For general r, one proceeds analogously and finds 

j fc+^-lAT'-l 

(15) 77((fc+^)m + r) = ^lim^ - £ £ ^ Mn+ ^ M(n+m)+r+t . 

^ ' t=0 n=0 

One now has to split the sum over t according to the ten cases of Table [TJ In each row, the 
condition for t is formulated in such a way that the difference of the bounds gives the proper 
multiplicity of the resulting term to the sum, which may be zero in some cases. 

Observing v n = — v n , one simply has to add the terms of the form v n v n+m with their signed 
multiplicities, which contribute to rj(m), and those of the form v n v n+m+ i, which contribute 
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to 77(777.+ 1). For instance, when < r < k, one finds the multiplicity of v n v n+m as 

(k—r) — (min(/c, k+£— r) — [k—r)j + (max(£;, k+i—r) — £;) 

= k + £ — r — 2(min(fc, k+£ — r) — /c + r) = k + £ — r — 2 min(r, £) 
= k + £ — r — 2 min(/c, £, r, k+£— r) = a klr 

where we used min(a, 6) + max(a, b) = a + b and, in the last line, the inequality < r < k. The 
required denominator (k + €) in the claimed recursion emerges from the splitting as shown 
above in Eq. (|15p . Likewise, the multiplicity for v n v n + m +i calculates as 

{k — min(/c, k+£— r)) — [(k+l) — max(/c, k+£— r)) 
= — £ + max(/c, k+£— r) — min(A;, k+£— r) 
= 2k — r — 2 min(fc, k+£ — r) = r — 2min(r, £) = a kek+i _ r , 

where we used that min(r, £) = min(fc, £, r, k+£ — r) holds in this case. The remaining cases 
(for k < r < k+£) follow from similar calculations; see Tabled! This completes the argument 
for 77(777) with 777 > 0. 

Finally, we know that r/(—n) = 77(77) for all 77 € 7L. Let m < be fixed, so that m = —\m\ 
with |?77| > 0. If r = 0, one simply has 

r][{k+£)m) = 77((/c+£)|777|) = 77(|?77|) = 77(777). 

When 1 < r < k+£, set m' = \m\ — 1 and s = k + £ — r, so that m! > and 1 < s < k+£. 
Then, one finds 

r]({k+£)m + r) = 77((A;+^)|t77| — r) = rj((k+£)m' + s) 

= ^-j( a kAsV( m ') + UkAk+e-sV(™' + l)) = j^rj (a kt e, r + a k,e,k+e-r v(m+l)) 

due to the reflection symmetry of 77 together with the recursion for positive arguments. □ 

Note that the recursion in Lemma [T] uniquely determines all coefficients 77(777) once 77(0) is 
given. Moreover, the recursion is linear, which will have strong consequences later. 

Since, for any given k,£ € N, every member of the corresponding hull of weighted Dirac 
combs has the same autocorrelation measure, we speak, from now on, of the autocorrelation 
measure of the gTM system (for parameters k,£ E N). Let us define £(N) = Yln=-N{ r l(. n )) > 
where we suppress the parameters k and £. For k = £ = 1, we know a bound on the growth 
rate of U(N), namely S(4iV) < |S(27V), from [3H For k + £ > 2, we formulate a similar 
result (with a technically more involved but structurally slightly simpler proof) as follows. 

Lemma 2. Let k,£ G N with k+£ > 2 be fixed, and let 77(77) with n € Z be the corresponding 
autocorrelation coefficients from Lemma [H Then, there is some positive number q < k + £ 
such that U((k+£)N) < qS{N) for all N € N. 
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Proof. Using the recursions of Lemma HJ one finds 

(k+e)N k+e-i N-i 

E((k+l)N)= ^{n)f = ( V ((k+£)N)) 2 + E mk+Dm + r))' 

n=-(k+£)N r=0 m=~N 

k+l-1 N-l 

= E(N) + + Y E ("M^M + a kAk+e _ r 7](m+l)) 2 



r=l m=—N 

E(N) ( 



S ( N ) (n „s2 2 2 A A 

i k + «) + 2^ ( a M,r + a k,e,k+e-r) j + 

, — 1 / 



(k + V ^ " r **" +t - r > ; 1 (k + £) 



r=l 

with ^ = (Em=-jv 2 k(m)r ? (m+l)|) (Erii~Vfe,£,r a k,e,k+e-r\) bem S a sum of non-negative 
terms only. Noting that 

N-l N-l 

Y 2\r](m)r](m + l)\ < ^ (?](m)) 2 + (^(m+l)) 2 < 2U(N), 

m=—N m=—N 

one obtains A < E(N) X^rii" 1 2 l a fc ir a k i k+£-r\- Employing the binomial formula results in 

z(( k +£) N ) < (|^((^+^) 2 + fc EVM,ri + i« fciW -,i) 2 )- 

Our claim follows if we show that the term in the large brackets is smaller than (A; + £)^ . For 
1 < r < k + 1 — 1, we know that 1 < min(/c, r, /c + ^ — r) < min(r, /c + £ — r), which implies 
\o k e r \ < k + £ — r and hence 

( 16 ) \ a k,£,r\ + \ a k,e,k+t-r\ <k + t. 

Since A; + £ > 2 by assumption, the stronger inequality \a k e 1 | + \a k e k+i _i\ < k + £ — 2 holds 
for r = 1, so that at least one term in the sum is smaller than k + £. This means that a 
q < k+£ exists such that E((k+£)N) < qE(N) holds for all N > 1. □ 

The recursion derived in Lemma Q] can now be used to show the absence of pure point 
components (by Wiener's criterion, which will rely on Lemma [2]) as well as that of absolutely 
continuous components (by the Riemann-Lebesgue lemma, which will rely on the special 
relation 77 ( ( A; + £)mj = r/(m) from Lemma [I]), thus establishing that each sequence in this 
family leads to a signed Dirac comb with purely singular continuous diffraction. 

Theorem 1. Let k,£ £ N. The diffraction measure of the gTM substitution g = g kl is 
the Fourier transform 7 of the autocorrelation measure 7 of Lemma [H It is the diffraction 
measure of every element of the hull of weighted Dirac combs for g. Moreover, 7 is purely 
singular continuous. 

Proof. Since the statement is clear for k = £ = 1 from Section [2] together with [32j H] , let 
k + £ > 2 be fixed. The corresponding autocorrelation is unique by Lemma HJ Since it is 
positive definite by construction, its Fourier transform exists [21], and then applies to each 
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element of the hull again. Since the underlying Dirac comb is supported on Z, we know from 
[U Thm. 1] that 7 is Z-periodic, hence it can be written as 

(17) ^y = fl *5 z with M = 7j [01) . 

Here, /j is a positive measure on the unit interval (which is a representation of the unit 
circle here), so that the inverse Fourier transform /x, by the Herglotz-Bochner theorem, is 
a (continuous) positive definite function on Z (viewed as the dual group of the unit circle). 
Since 7 = ft 5% by the convolution theorem together with the Poisson summation formula 
5z = $z, we see that this function is 

(18) 77(777) = f e 2nimx dv(x) = /2(m). 

Jo 

Let us now use the recursion for 7/ to infer the spectral nature of [/, and thus of 7. 

Lemma [2] implies -^S(JV) N ^°°) 0, which means (7) pp = by Wiener's criterion [53] ; see 
also [33l|43]. We thus know that 7 = (7) + (7) is a continuous measure, where the right- 
hand side is the sum of two positive measures that are mutually orthogonal (in the sense that 
they are concentrated on disjoint sets). Each is the Fourier transform of a positive definite 
measure with support Z, hence specified by autocorrelation coefficients rj sc and r/ 3C which 
clearly satisfy rj sc (m) + rj ac {m) = 77(771) for all m £ Z. The recursion relations for 77 imply a 
corresponding set of functional relations for the non-decreasing and continuous distribution 
function F defined by F(x) = 7Q0, x]) for < x < 1. Due to the orthogonality mentioned 
above, the same relations have to be satisfied by the ac and sc parts separately. This in turn 
implies that 77 sc and ?7 ac must both satisfy the recursion relations of Lemma [TJ however with 
a yet undetermined value of 7? ac (0), and ?7 SC (0) = 1 — i] ac (0). 

The recursion of Lemma Q] with m = and r = 1 gives 

while r = leads to 77 ac ((/c+^)m) = f] ac (m) for all m 6 Z. Since we have lim^^oo i] ac (n) = 
from the Riemann-Lebesgue lemma, compare [33J, we must have 7/ ac (m) = for all m > 0. 
When k + £ > 3, 77 ac (l) = forces 77 ac (0) = 0, and then r] ac (m) = for all m G Z by the 
recursion, which means (7) = 0. When k + t = 3, we have 7/ ac (l) = 0, but the recursion 
relation for m = and r = 2 leads to 77 ac (2) = — | ?7 ac (0), hence again to r] ac {0) = with the 
same conclusion. 

As a consequence, ?7 SC (0) = 1 and i] = rj sc . We thus have 7 = (7) as claimed. □ 

Remark 2 (Diffraction with general weights). If an arbitrary gTM sequence is given, the 
diffraction of the associated Dirac comb with general (complex) weights h± can be calculated 
as follows. If h is the function defined by h(l) = h + and h(l) = h_, one has 

E11 \ c h + +h_ h + -h_ 
h(w n )5 n = 6 Z + u 
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with the oj from Equation (|13p . The autocorrelation of u, clearly exists and calculates as 

\h,+h_\ 2 \h,-h_\ 2 
Ih = 4 h + 4 7 

with 7 from (|14p . This follows from 5% = 5 Z together with <5 Z © oj = 5 Z ® oj = 0, which 
is a consequence of the fact that 1 and I are equally frequent in all gTM sequences. The 
diffraction is now obtained as 

_ |/i + +/i_| 2 e \h + -h_\ 2 ^ 
In = 1 <*z + Z 7 

by an application of the Poisson summation formula 5 Z = <5 Z . Since 7 is purely singular 
continuous, this is a diffraction measure with singular spectrum of mixed type. 

This diffraction does not display the full dynamical spectrum of the gTM system, which is a 
well-known phenomenon from the classic TM system |25j . In the latter case, this is 'rectified' 
by the period doubling system as a topological factor. We will return to this question for the 
gTM systems in Remark HI 



4. The diffraction measure of the gTM system 

Let us consider the diffraction measure in more detail, which we do via a suitable dis- 
tribution function F for the (continuous) measure 7. This is done as follows. First, we 
define F{x) = 7QO, 1]) = /u([0, x]) for x G [0,1]. The Z-periodicity of 7 together with 
/u([0, 1]) = 77(0) = 1 means that F extends to the entire real line via F(x + 1) = 1 + F(x). 
Moreover, since 7 = 7, we know that 7 is reflection symmetric. With -F(O) = 0, this implies 
F(—x) = —F(x) on R, which is our specification of F in this case. 

Proposition 2. Let k,£ £ N be fixed. The distribution function F of the corresponding 
diffraction measure is non- decreasing, continuous, skew- symmetric and satisfies the relation 
-F(x + 1) = 1 + -F(x) on the real line. Moreover, it possesses the series expansion 

\ - rj(m) 

Fix) = x + > sm(27rmx), 

^ run 

m>l 

which converges uniformly on M. 

Proof. By construction, F is non-decreasing, and is continuous by Theorem [TJ So, F(x) — x 
defines a 1-periodic continuous function that is skew-symmetric and the difference of two 
continuous, non-decreasing functions, hence it is of bounded variation. By standard results, 
see [441 Cor. 1.4.43], it has thus a uniformly converging Fourier series expansion 

00 

F(x) — x = b m sin(27rmx). 

m=l 
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The Fourier coefficient b m (for m G N) is 

b m = 2 I sm(2irmx) (F(x) — x) dx = + 2 / sin(27rmx) F(x) dx 

Jo rn-ir J 

= — f cos(2nmx) dF(x) = — C e 2 ~ dF(x) = 

The first step in the second line follows from integration by parts, while the next is a con- 
sequence of the symmetry of dF together with its periodicity (wherefore the imaginary part 
of the integral vanishes). Recalling that F (restricted to [0, 1]) is the distribution function of 
the probability measure fi completes the argument. □ 

It is interesting that the autocorrelation coefficients occur as Fourier coefficients this way. 
In preparation of a later result, let us look at this connection more closely. A key observation 
is that the recursion relations for n from Lemma [T] (which have a unique solution once the 
initial condition 77(0) = 1 is given, with \r](n)\ < 1 for all n G Z) can also be read as a recursion 
as follows. Let fi G [—1, 1] N be a sequence and define a mapping \E' via 

a k,£,r + a k,e,k+e-r Pn+i > if 7i = and 1 < r < k + £, 
( 19 ) ^P) {k+ l)n + r ■= Fb \ {k + l) Pn ' if n G N and r = 0, 

. a k,£,r Pn + a k,i,k+£-r Pn+i , if re G N and I < r < k + £, 
which completely defines the sequence *$>(3. This mapping derives from the recursion for 77 
with positive arguments when r/(0) = 1. 

Lemma 3. The mapping ^ maps [— 1,1] into itself, with W^SPW^ < Wf^W^. Moreover, 
for any G [-1,1] N , the iteration sequence pW defined by ^ N+1 ^ = for N > 

converges pointwise towards (?7( n )) ne ^- 

Proof. The first claim follows from |a fc ^ r | + \a k e k+e _ r \ < k + £, which was used earlier in 
the proof of Lemma [2] via the triangle inequality. When r = or when k = I = r, one has 
equality here, so that \E f is not a contraction on [—1, 1] N for the supremum norm. 
Observe that the iteration for f3 t is closed and reads 

a, _ k + e-3 1 o 

Pi ~ k + £ — £ Pi ' 

which is an affine mapping with Lipschitz constant and hence a contraction. The iteration 
for /3 1 thus converges exponentially fast (to fxfef ) by Banach's contraction principle. 

What happens with the iteration for /3 1 determines everything else, because the components 
fi n with (k + £) m < n < (k + £) m+1 and m > emerge from fi 1 in m steps of the iteration. 
In particular, the iteration also closes on any finite block with 1 < n < (k + £) m and fixed 
m G N, and shows exponentially fast convergence. Note though that the iteration is only 
non-expanding as soon as m > 1, while \& induces an affine mapping with Lipschitz constant 

_ max{\a kAk+e _ r \ \ 1 < r < k + £} k + £-l 
L ~ kTl " k + £ < lj 
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on the components (3 n with 1 < n < k + 1. 

Pointwise convergence is now clear, and the limit is the one specified by the original recur- 
sion, which proves the claim. □ 

The recursion relations for 77 can also be used to derive a functional equation for the 
distribution function F. Observe first that 

rl rk+£ 

r]((k+£)m + r) = / e 2ni « k+ V m+r > dF{x) = e 2nimx e 2 ™^ dF(^) 
Jo Jo 

( 2 °) , 

= / e 2 ~ e 27ri ^R e 2 ™^dF(f±§). 

J° s=0 

On the other hand, we know from Lemma [T] that 

rj((k+£)m + r) = -^-j (a k/r r?(m) + a k/jk+i _ r rj(m+l)) 

(21) n ^ 2nix 

2mmx a k,i,r + a k,e,k+£-r e , „, , 

o 6 k + e a (x) - 

A comparison of (|2Up and (|21|) leads to the following result. 

Proposition 3. The distribution function F for k, £ G N satisfies the functional equation 
Ak+t)x 2 k+e-l 

= ~kTI I ^(^) dF ^ with = 1 + ~k~TI 2 ° k ^ r cos ( 27Trx ) ■ 

r=l 

This relation holds for all x E R, and 1? is continuous and non-negative. 



Proof. Equations (J20J) and (|2ip . which hold for all m € Z, state the equality of the Fourier 
coefficients of two 1-periodic Riemann-Stieltjes measures, which must thus be equal (as mea- 
sures). For < r < k+l, we thus have 

s=0 

Fix an integer i with < i < A; + £ and multiply the equation for r on both sides by 
exp(— 27ri-j^) ■ Since ^r=o -1 ex P (^ 7r ^ r fc+i ^ ) = (k+-0^a,t> a summation over r followed by a 
division by (fc + leads to 



which is valid for all x £ [0, 1). 
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To derive the functional equation, we need to calculate F(x) and hence to integrate the 
previous relations with an appropriate splitting of the integration region. When [y] and {y} 
denote the integer and the fractional part of y, one finds 

(22) F(x) = / dF(y) = / dF(^M) + £ / dF(gf) , 

J o J o t=0 Jo 

which holds for all x G [0, 1). Observe next that 

d^(iS) = kTil' H^)^(y) = ^^%fe)dF(z). 

holds for any < t < k + £, where we used that dF(z — n) = dF(z) for all n G Z as a 
consequence of the relation F(z + 1) = 1 + F{z) for z£l. Similarly, one obtains 

r {(k+l)x} .. r (k+£)x 

/ dF(^M)= 1 / *G&)dF(,). 

We can now put the pieces in (|22p together to obtain the functional equation as claimed, 
which clearly holds for all x G R. 

The continuity of $ is clear. Its non-negativity follows from the functional equation, because 
we know that F is non-decreasing on [0, 1] (as it is the distribution function of the positive 
measure [/,). If we had i?(a) < for some a G [0, 1], there would be some e > such that 
"&{y) < —e in a neighbourhood of a, which would produce a contradiction to the monotonicity 
of F via the functional equation. □ 

Remark 3 (Properties of the integration kernel). The non- negative function •& of Propo- 
sition [3] has various interesting and useful properties. Among them are the normalisation 
relations 



f 1 1 

/ $(x)dx = 1 and / fi(x)dx = — , 

Jo Jo 2 



which follow from Jq cos(27rrx) dx = for r ^ together with the 1-periodicity of i? and its 
symmetry (whence we also have — x) = "!?(x)). Another is the bound 

Moo < ?, 

where g is the number from Lemma [2] (when A; + ^ , >2)orq , = 2 (when fc = £ = 1). This 
bound is proved by another use of Eq. (fT6j) and the lines following it. 



The functional equation of Proposition [3] provides the basis for the calculation of F by a 
Volterra-type iteration. To this end, one defines F (x) = x (so that dF (x) = dx) together 
with the recursion 
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for n > 0. It is clear that each F n defines an absolutely continuous Riemann-Stieltjes measure, 
so that one can define densities /„ via dF n (x) = f n (x) dx. This gives 



f n+1 (z)dz = F n+1 (x) = / 0(z)f n {(k + e)z)dz, 

Jo 

which results in f n+ \(z) = $(z) f n {(k + £)z), and hence in the continuous function 

n-1 

(24) f n (z) = l[V((k + £) j z). 

j=0 

To put this iteration into perspective, let us introduce the space D of non-decreasing and 
continuous real- valued functions GonR that satisfy G(—x) = —G(x) and G(x+1) = 1 + G[x) 
for all x € R. In particular, this implies G(q) = q for all q € ^Z. We equip this space with 
the 11.11 oo-norm, and thus with the topology induced by uniform convergence, in which the 
space is closed and complete. Each G G D defines a positive Riemann-Stieltjes measure on 
K that is reflection symmetric and 1-periodic. Also, G(x) — x always defines a continuous, 
skew-symmetric and 1-periodic function of bounded variation. Our distribution functions F 
from above are elements of D. 

Let k, £ £ N be fixed. Define a mapping $> by G H > $G, where 

r(k+e)x 

(25) = -L_ ^ H^- e )dG(y). 

Clearly, our previous iteration (|23p can now be written as F n+ \ = $>F n with the initial 
condition F Q (x) = x, where F € D. 

Lemma 4. The operator $ maps D into itself. Moreover, for arbitrary F^ € D, the iteration 
sequence defined by F^ N+l "> := $F( N ' for N > converges uniformly to the continuous 
distribution function F of Proposition [2l 

Proof. Let G G D. It is clear that $G is again continuous (since i? and G are continuous) 
and non-decreasing (since •& is non- negative by Proposition [3]). The skew-symmetry follows 
from dG(— y) = —dG(y) by a simple calculation. Finally, one has 

r (k+e)(x+i) 

$ G{x + i) = _±_ J (_IL_) dG{y) = $ G (x) + / , 

where the remaining integral /, using $ from Proposition [3] and the periodicity of dG, is 



fc+^-1 f t+i „ r l 



fc+«-l „ t+ i ,1 



t=0 " 1 t=0 
jfe+i-1 r l k+i-r 



1 E «M,r / E cos(2vr^)dG( y ) = 1. 

r=l 1/0 t=0 



_2 
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The last step follows because the sum under the integral vanishes as a consequence of the 
relation YlT=o ex P 

(2m £f) = for 1 < r < m. 
To establish the convergence, let us first show pointwise convergence 
for all x G R. Observe that F^ N \x) — x defines a 1-periodic, continuous function on R for each 
N G No, and so does F(x) — x. Their restrictions to [0, 1] define regular (signed) measures on 
the unit circle, i/W and v say, with Fourier-Stieltjes coefficients a n N ^ and a n , where n G Z. 
The latter coefficients follow from Proposition [21 and the former from the observation that 
every G G D has a uniformly converging Fourier series of the form 

G(x) = x + — sin(27rnx) , 
^-^ nir 

n=l 

where \/3 n \ < 1 by an application of [56, Thm. II. 4. 12]. Here, we have /3( w+1 ) = with 
the mapping \& of Lemma [3l with initial condition (3^ and convergence /% jV ~ ! ' 00 > r](n) for 
each n G N. Consequently, a n N ~ Jf00 > a n for all n G Z, which means that ^W(p) ^^(p) 
for any trigonometric polynomial p with period 1, and hence (by an application of the Stone- 
Weierstrass theorem [46J) for all continuous functions on [0, 1]. This proves weak convergence 
V {N) Af ~ > °° > v on unit circle. Since all measures are absolutely continuous, with continuous 
Radon-Nikodym densities, this implies pointwise convergence F^ N \x) — x N ^°°) F(x) — x 
for all x € [0, 1] , and hence (by periodicity) the pointwise convergence on R claimed above. 

Uniform convergence on R now follows from that on [0,1], which can be shown via the 
'stepping-stone' argument from the proof of [19^ Thm. 30.13]. For completeness, we spell out 
the details. Given e > 0, there are numbers m G N and = x < x 1 < ■ ■ ■ < x m = 1 with 

\F( Xi ) - F(x^t)\ = F( Xi ) -F(xi-i) < e 

for 1 < i < m. Also, for all sufficiently large N € N, one has 

\F^( Xi )-F( Xi )\ < e 

for all < % < m. Since = F(l) = 1, consider now an arbitrary x £ [0,1), so 

that Xi-i < x < Xi for precisely one i G {1,2, . . . ,m}. Using monotonicity, this implies the 
inequalities 

F(x<_i) < F(x) < F{ Xi ) < F(xi_i) + £ 

and 

F(x<_i)-e < fW(^i) < FW(i) < FW( X .) < F(x 4 )+e < F( Xi -i) + 2e. 

Together, they give |F^(x) — F(x)| < 2e, which holds for all x G [0, 1], and then uniformly 
for all x G R, as - F is 1-periodic for all N G No- □ 

Due to our convergence results, Equation (|24j) means that the measure 7 has a (vaguely 
convergent) representation as the infinite Riesz product Y[ n >o "&{{k+() n x) . The entire analysis 
is thus completely analogous to that of the original TM sequence and shows that the latter 
is a typical example in an infinite family. Two further cases are illustrated in Figure [3l 



20 MICHAEL BAAKE, FRANZ GAHLER, AND UWE GRIMM 




Figure 3. The continuous and strictly increasing distribution functions of the 
generalised Thue-Morse measures on [0, 1] for parameters (k, £) = (2, 1) (left) 
and (5, 1) (right). 

Remark 4 (Pure point factors). The block map ([8]) applies to any member of the gTM 
family, and always gives a 2-to-l cover of the hull X^ of the generalised period doubling 
(gpd) substitution 

(26) Q=Qk,l- a^b^ab^b, b ^ b^ab^a. 

Since we always have a coincidence (at the kth. position) in the sense of Dekking [24], they all 
define systems with pure point spectrum (which can be described as model sets in the spirit 
of [16\ I15j ) - another analogy to the classic case k = £ = 1. Also, for given k,£ G N, the 
dynamical system (X^,Z) is a topological factor of (X™,Z). 

Moreover, the dynamical spectrum of the gTM system contains Z[^jy as its pure point 
part, which happens to be the entire spectrum of the gpd system. The latter fact can be 
derived from the support of the gpd diffraction measure, via the general correspondence 
between the dynamical and the diffraction spectrum for pure point systems [371 [12] . The 
detailed calculations can be done in analogy to the treatment of the period doubling system 
in [151 E], after an explicit formulation of the one-sided fixed point of the gpd substitution, 
which results in the Fourier module Z[^|]. Consequently, the gpd system is a topological 
factor of the gTM system with maximal pure point spectrum. 

The gpd system can be described as a model set (with suitable p-adic type internal space). 
From [13], we know that there exists an almost everywhere 1-to-l 'torus' parametrisation 
via a solenoid Sk+e- Here, the inflation acts as multiplication by k+i. In fact, the solenoid 
provides the maximum equicontinuous (or Kronecker) factor of the gTM system. 



GENERALISED THUE-MORSE SEQUENCES 



21 



5. Topological invariants 

The formulation via Dirac combs embeds the symbolic sequences into the class of translation 
bounded measures on R, as described in a general setting in [12]. It is thus natural to (also) 
consider the continuous counterpart of the discrete hull X in the form 

Y = {S t *u | t € R}, 

where u is an arbitrary element of X (which is always minimal in our situation) and the 
closure is now taken in the vague topology. Here, Y is compact, and (Y, R) is a topological 
dynamical system. Note that the continuous tiling hull mentioned earlier is topologically 
conjugate, wherefore we use the same symbol for both versions. The discrete hull X, in this 
formulation, is homeomorphic with the compact set 

Y = {v G Y | K{0}) + 0} C Y, 

while (Y,R) is the suspension of the system (Y ,Z). 

We are now going to construct the continuous hulls of the generalised Thue-Morse and 
period doubling sequences as inverse limits of the inflation map acting on an AP-complex 
[3], and use this construction to compute their Cech cohomology. At first sight, there are 
infinitely many cases to be considered. We note, however, that the AP complex T depends 
only on the atlas of all r-patterns occurring in the tiling, for some bounded r. In the original 
AP construction [3], T consists of 3-patterns - tiles with one collar tile on the left and one 
on the right. How these collared tiles are glued together in the complex is then completely 
determined by the set of 4-patterns occurring in the tiling. Therefore, there are only finitely 
many different AP complexes T to be considered. 

For obvious reasons, we want to extend the factor map <p : Y™ ->■ Y p k d e between the 
continuous hulls to the cell complexes approximating them, compare (|1Q[) . We therefore 
choose the cell complex T™ to consist of tiles with one extra layer of collar on the right, 
compared to the cell complex r?^, so that we obtain, for each pair (k,£), a well-defined factor 
map (/>: T™ — > r?*^, which we denote by the same symbol. 

While the extra collar seems to complicate things, we can compensate this by a simplifi- 
cation compared to the original AP construction. In [29], it was shown that a computation 
using a complex with one-sided collars already yields the correct cohomology groups, even 
though the inverse limit is not necessarily homeomorphic to the continuous hull of the tiling. 
The minimal setup therefore consists of AP complexes T^^ 1 with tiles having two-sided col- 
lars, and complexes r?^ having tiles with left collars only. All collars have thickness one. A 
straight-forward analysis shows that only three cases have to be distinguished: k = £ = 1 (the 
classical case), either k = 1 or £ = 1 (but not both), and k, £ > 2. The complexes T™ are 
all sub-complexes of the complex shown in Figured! If k = £ = 1, the two loops on the left 
and right have to be omitted, because the corresponding patterns 111 and 111 do not occur 
in the classic TM sequence. Likewise, in the case k,£ > 2, the lens in the centre has to be 
omitted, whereas the full complex has to be used in all remaining cases. 
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Figure 4. The universal complex for the gTM system. Edges are labelled 
both with the corresponding collared tile of the gTM tiling and its image in 
the gpd tiling. 



ba 




ab 



Figure 5. The universal complex for the gpd system. Its edges correspond 
to tiles with one-sided (left) collar. 

For the gpd sequences, an approximant complex that includes all different cases as sub- 
complexes can be constructed in a similar way. This complex is shown in Figure [5j As we 
have labelled the edges in Figure |4] with both the 3-pattern in the gTM sequence and the 
corresponding image under (ft in the gpd sequence, it is easy to see how the factor map (ft maps 
the complex of Figure [4] to that of Figure [5j the loops on the left and right of Figure |4] are 
both mapped onto the right loop of Figure the central lens in Figure [4] is wrapped twice 
around the left loop in Figure \5\ and the rhombus in Figure S] is wrapped twice around the 
central loop in Figure Therefore, consists of the left and the central loop of Figure [5] 
in the case k = £ = 1 (the classical case), the central and the right loop of Figure [5] in the 
case k,£ > 2, and of all three loops in the remaining cases. It is obvious that the map (ft is 
uniformly 2-to-l also at the level of the approximant complexes. 

The cohomology of the continuous hulls and is now given by the direct limit 

of the induced maps p* kl and p ki on the cohomology of the AP complexes, i?*(r™) and 
H*(T^). While there are only 3 different cell complexes in each case, there is an infinite 
family of maps p* k t , respectively p k£ , but these can be parametrised by k and t. 
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We first define a basis of the homology of the full complexes of Figures H] and subsets 
of which shall be used for all k and £, and pass later to the corresponding dual basis for the 

TM 

kl 



cohomology. The homology of the full complex T T ¥ is generated by the basis 



c™ = m, c™ = m, c™ = im, c™ = i, c™ = i, 

where the words on the right have to be thought of as being repeated indefinitely. Likewise, 
the homology of r£ d is generated by the cycles 

c^ d = a, c?. d = ab, Cg d = b. 

Clearly, we have 

Uc? M ) = <? + cf, McI M ) = cf + 4 d , McJ M ) = 2cf, 

McJ M ) = cf , McJ M ) = cf. 

Here, the lower asterisk in denotes the induced action on homology, whereas (jf denotes 
the action on cohomology. 

In order to go into more detail, we have to distinguish the three different cases. We begin 
with the classical case, k = £ = 1. Here, the relevant basis elements are c TM , c™ and c™ 
for TM, and c pd and c^ d for pd. On these, the substitution acts as 

P*(C™) = C™ + C™, p,(c™) = C™ + C™, p,( C ™) = 2(C™ + C TM) _ C T M) 

p^cf) = cf, p'Acf) = 2cf + cf, 

where and p'^ again denote the induced action on homology. If we express these maps 
as matrices A™ and A pd acting from the left on column vectors, with respect to the basis 
above, and likewise define a matrix P for the action of we obtain 



(27) A pd = , A 




iTM 




A i 

11 2 | . P 

^0 

The corresponding action on the dual basis of cohomology is simply given by the transposed 
matrices, or by the same matrices regarded as acting from the right on row vectors. We adopt 
the latter viewpoint here. It is easy to verify that the matrices (|27p satisfy A pd P = P^4™. 
In other words, the substitution action commutes with that of A pd and ^4™ have left 
eigenvectors and eigenvalues 

(28) _ 2 x j and 

Next, we look at the case k,£ > 2, which is still relatively simple. Here, we have to work 
with the basis elements c™, cJ M and c™ for gTM, and c^ d and Cg d for gpd. The action of 
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the substitution is then given by 



4 J 
TM\ 



P*(cJ M 
P*(cJ 

p'Acf) 

p'M d ) 



3c™ + (2(k + e)-(5)(c™ + c™), 

2)c™, 



TM 



3 C P d + 2{k + i-3)cf, 



TM 



2cf + {k + £- 4)d 



Pd 
3 • 



The corresponding matrices , ^4™ and P k e read 



A pd _ 3 2 \ 4 tm 

\2(k+£-3) k+i-Aj ' 



3 1 1 \ 

2(k+£)-6 k-2 1-2 



p -( 2 









1 





and satisfy A^ d e P k e = P kl A™. The matrices A k d e and A™ have left eigenvectors and 
eigenvalues 



(29) 





and 




1 1 | k+l\ 
-1 -1 | -1 
1 -1 | k-i / 



Finally, the case where min(fc,£) = 1, but k + £ > 2, requires the full AP complex. For 
notational ease, we consider the case k = 1, i > 2. The case i = 1, A; > 2 is completely 
analogous. The substitution then acts as follows on our basis: 



P*(c? M ) 


= C JM + C TM + 1)c TM 


+ {21- 


3)cJ M 


P*(cJ M ) 


= C ™ + c TM + (£ _ 1)c TM 


+ {2£- 


3)c™, 


P*(cJ M ) 


= c TM + c TM + c TM + (2 ^. 


-3)(c™ + C ™) 


P*(cJ M ) 


= c™ + (£-2)c™ 






P*(cJ M ) 


= c™ + (f-2) C ™, 






p'M d ) 


= c f + (£-l)cf, 






p'M d ) 


= cf + 2cf + {2£-3)cf, 






p'M d ) 


= cf + cf + {£-2)cf. 
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2t-3 1-2) 



1 



1 



The corresponding matrices are 

t ° 

/ 1 1 \ 1 

1 2 1 , A™= 1 

i-\ 21-3 21-3 
\2£-3 £-1 2£-3 

/l 1 
112 
^0 1 1 

In this case, A™ and have left eigenvectors and eigenvalues 

/ 3 3 4 1 




1 




£-2 



1 \ 




-2 
/ 



1 



(30) 



1 1- 




and 



1- i l-£ 
1 -1 

2- 1 2-1 



2-1 1 

1 

2-2£ 1 

1 



1 
1 

-1 
1 

-1 



1 



+ 1\ 
-1 



o / 



\2-£ 1-2 

As the cohomology of the continuous hull is given by the direct limit of the action of the 
substitution on the cohomology of the AP approximant complexes, it is determined by the 
non-zero eigenvalues of the matrices A™ and A^ d t For all k and £, these eigenvalues are k + £ 
and —1 for the gpd case, and k + £, —1 and k — I in the gTM case, where the last eigenvalue 
is relevant only if k ^ £. Since H°(T) = Z if T is connected, and the substitution action on 
H (r) is trivial, with eigenvalue 1, we can summarise our cohomology results as follows. 

Theorem 2. The Cech cohomology of the continuous hull of the gTM sequences is given by 
H°{Y™) = Z and 



z^jezez^], if\k-e\>2, 
z[^]ez 2 , if\k-£\ = i, 

z[^]e2, ifk = e. 



The Cech cohomology of the continuous hull of the generalised period doubling sequences is 
given by 

H°(Yg) = Z, ff 1 (Yg) = Z[ ] ^]©Z > 
valid for any pair k, £ G N. □ 

Since the cohomology of the hulls is given by the direct limit of the columns of the diagram 

H k Ff}) ^(rg) # fe (r™) 

(31) x(fc+£) e « e * 



+ ^ fe (r fc 7) 



2(i 
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(compare (|10p). we now have access also to the homomorphisms in the sequence 

(32) H*(S k+e ) H*(Y P k d e ) H*(Y™). 

For all three spaces, H° = Z, and the maps (ff and ip* acting on them are isomorphisms. 
Further, it is easy to see that tp* embeds .ff 1 = Zf^^j] non-divisibly in H l (Y 1 ^ d i ); it is 

simply mapped isomorphically onto the summand Z[^|] in ff 1 (Y^). The same cannot be 
said of the second map, however. For all pairs k, £, the matrix P maps the eigenvector of A? . 
with eigenvalue —1 to twice the corresponding eigenvector of ^4™. As a result, the quotient 
Hl ( Y ™) /p*(H l (Yl d e )) has a torsion component Z 2 ; the summand Z of H l (Y p k d ) comine from 
the eigenvalue —1 is mapped to 2Z in if 1 (Y™). This result was already known for the classic 
period doubling and Thue-Morse sequences |18j . and extends to the generalised sequences. 
Therefore, in the case k = £, where the cohomologies of gTM and gpd are isomorphic as 
groups, one should rather regard i? 1 (Y^) as subgroup of index 2 of if 1 (Y™). 

6. Dynamical zeta functions 

The continuous hull permits to employ the Anderson-Putnam method [3] for the calculation 
of the dynamical zeta function of the inflation action on Y. The dynamical zeta function [IS] 
of a substitution can be viewed as a generating function for the number of fixed points a{n) 
under n-fold substitution via 

OO , s 

(33) C(,) = exp(£^4 

n=l 

Knowing the fixed point counts a(n), one can calculate the number c(n) of cycles of length n 
from the formula 

1 



d\n 

where /i is the Mobius function from elementary number theory (and should not be confused 
with the measure \i that appeared above). 

Anderson and Putnam [3] showed how the dynamical zeta function of a substitution tiling 
can be expressed by the action of the substitution on the cohomology of the AP complex T. 
In the one-dimensional case, the zeta function is given by 

det(l - zA°) 
C(Z) ~ det(l - zA 1 ) ' 

where A k is the matrix of the substitution action on H k (T). In our case, A is a 1 x 1 unit 
matrix, and A 1 is diagonalisable, so that we can rewrite the zeta function as 

C{z) = ILU-W 

where Aj are the eigenvalues of A 1 . For the gTM and gpd tilings, these eigenvalues have been 
derived above, so that we arrive (after a simple calculation) at the following theorem. 
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Theorem 3. Let k,£ € N. The generalised Thue-Morse sequence defined by the inflation g ki 
°f (llip possesses the dynamical zeta function 

tilt*) ' ~ " 



{i + z ){i-{k + i)z)(i-{k-i) z y 

while the induced generalised period doubling sequence, as defined by g' ke of ([26]) . leads to 



(l + z)(l - (k + £)z) 

In particular, when k = £, one has £™ = £ pd . □ 

For our two systems, the corresponding fixed point counts, for n E N, can now be calculated 
from Eq. (f33l) to be 

(34) af/n) = (k + If - (l - (-1)") and a™( n ) = ag(n) + (fc - <)» 

It is an interesting exercise to relate the orbits of the two systems according to the action of 
the mapping in agreement with these counts. 

There is an interesting connection between the inflation on or Y?^ and the multiplica- 
tion by (k + £) on the (matching) solenoid Sk+e, which emerges from the torus parametrisation 
|10| . 150 1 ^b3]- This solenoid is a set on which the multiplication is invertible. It is constructed 
via a suitable inverse limit structure (under iterated multiplication by (k + £)), starting from 
the 1-torus (or unit circle), represented by the unit interval [0, 1) with arithmetic taken mod 
1. Together with the above, we obtain the following commutative diagram 



(35) 



x(k+i) 



where ip denotes the torus parametrisation for the generalised period doubling sequence, in 
the spirit of [50^ [131 116j . The mapping ip is 1-to-l almost everywhere. Like for the classic 
period doubling sequence, it fails to be 1-to-l on exactly two translation orbits, which are 
mapped to a single orbit. 

Counting finite periodic orbits under the multiplication action on the solenoid, however, 
means that the inverse limit is not needed, so that the corresponding dynamical (or Artin- 
Mazur) zeta function coincides with that of the toral endomorphism represented by multipli- 
cation by m = k + £ > 2. This, in turn, is given by 



l — mz 

by an application of |11[ Thm. 1]. The number of fixed points is given by a^(n) = m T ' 
where m = k + £ as before for a comparison with (|3i 
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7. Further directions 

A natural question concerns the robustness of the singular continuous spectrum under 
simultaneous permutations of positions in g(l) and g(l), as considered in |55| . For k = I = 1, 
the two possible rules are the Thue-Morse rule (11, 11) and its partner (11, 11), written in 
obvious shorthand notation. Since the squares of these two rules are equal, they define the 
same hull, and hence the same autocorrelation. 

There are three possible rules for k = 2 and I = 1, namely (111, 111), (ill, III) and 
(111, 111). The first is our Q2 1 from above, while the third results in a fixed point that is 
the mirror image, and hence possesses the same autocorrelation. The middle rule, however, 
has the periodic fixed point ... 111111 .. . and hence autocorrelation 7 = (5q — 5±) * 62% ■ Its 
Fourier transform reads 

7 = S z+ i , 

which is a periodic pure point measure. By general arguments, one can see that the diffraction 
measure of the balanced situation (with 1 and 1 being equally frequent) must be of pure type, 
and cannot be a mixture. On the basis of the results from [26, 27J, it is then clear that, given a 
permutation, the diffraction is either pure point or purely singular continuous. It remains an 
interesting question to decide this explicitly for general k and £, and a general permutation. 

Moreover, it is clear that similar structures exist in higher dimensions. Indeed, starting from 
the treatment of bijective lattice substitution systems in [26l [27] , it is possible to demonstrate 
the singular continuous nature for bijective substitutions with trivial height lattice and a 
binary alphabet, and to calculate it explicitly in terms of Riesz products. Details will be 
explained in a forthcoming publication [9]. 
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